U(l) Gauge Theory of the Hubbard Model : Spin Liquid States and Possible 
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We formulate a U(l) gauge theory of the Hubbard model in the slave-rotor representation. From 
this formalism it is argued that spin liquid phases may exist near the Mott transition in the Hubbard 
model on triangular and honeycomb lattices at half filling. The organic compound k — (BEDT — 
TTF)iCui (CN)3 is a good candidate for the spin liquid state on a triangular lattice. We predict a 
highly unusual temperature dependence for the thermal conductivity of this material. 
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The Hubbard Hamiltonian is the simplest model which 
exhibits the Mott transition. For small Coulomb repul- 
sion the Fermi liquid phase occurs. In the large Coulomb 
repulsion limit the model reduces to the t-J model which 
has insulating ground state and spin order at half filling. 
In the intermediate region a spin liquid phase, i.e. an 
insulator without spin order, may arise. This spin liquid 
phase is more likely to occur if there is frustration. 

Recent experiments indicate that the spin liquid 
ground state may be realized in the organic compound 
k- (BEDT-TTF ) 2 Cu 2 (C N)^which is just on insulat- 
ing side of the Mott transition^. While the high tem- 
perature spin susceptibility can be fitted to a Heisenberg 
5 = 1/2 model with an exchange energy of 250K, the 
spins do not order down to 32mK. The susceptibility 
and Knight shift 2] drop sharply below ~ 20K, but re- 
mains finite down to the lowest temperature. The Car- 
bon nuclear spin relaxation rate 1/TjT also saturates to 
a finite value at low T. An interpretation in terms of a 
spin liquid on a frustrated lattice has been proposed^ 
but an alternative explanation in terms of localization 
has also been suggested 2]. This system is effectively 
described by a Hubbard model at half filling on a 2- 
dimensional triangular lattice. On this lattice the Heisen- 
berg model is known to have an antiferromagnetic (AF) 
ground state. Recently Motrunich has shown using pro- 
jected trial wavefunctions that a spin liquid state with 
spinon Fermi surface may be stable if the t-J model is 
extended to include higher order virtual hoppings, such 
as the ring exchange termsQ- If true, this will be a 
realization of the resonating valence bond (RVB) idea of 
AndersonQ. In order to include such charge fluctuations 
effects, it is clearly desirable to study the Hubbard model 
near the Mott transition. 

There are numerical evidences that the Hubbard model 
on the triangular lattice does not have AF order for 
moderate Coulomb repulsion in the insulating phase[j|. 
However, the nature of the disordered state is not un- 
derstood. Recently Florens and Georges^ introduced 
the slave-rotor representation of the Hubbard model. 
This method is more economical than the conventional 
representation which requires 4 slave bosons. From 
this theory the Mott transition on the square lattice was 



successfully described at the mean- field level gj. Thus it 
is of great interest to see whether the spin liquid phase 
is predicted by the slave-rotor theory in the triangular 
lattice. If a spin liquid mean-field state is found, the 
next question is the stability of the state. Although it 
is expected that the U(l) gauge theory will emerge as a 
low energy theory, it is not clear how to derive it from 
the mean-field theory of Ref. 6]. The first objective of 
the present paper is to formulate a U(l) gauge theory 
of the Hubbard model in the slave-rotor representation. 
Then we apply the formalism to the triangular and hon- 
eycomb lattices to find the spin liquid mean-field states. 
The low energy theories of some spin liquid states are 
shown to be a compact U(l) gauge theory coupled with 
gapless spinons. Finally we discuss the deconfmcmcnt 
phase of the U(l) gauge theory in connection with the 
experimentally observed spin liquid behavior in the or- 
ganic material. 

The Hubbard Hamiltonian reads 



H = 



<i,j>,& 



(1) 

Here Ci a is the annihilation operator of electron with spin 



a at site i. = \U 



-Mi 



is the complex hopping inte- 



gral with Aij representing either external electromagnetic 
(EM) vector potential or an intrinsic phase coming from 
the overlap of atomic wavefunctions. U is the on-site 
Coulomb repulsion. In the slave-rotor representation the 
electron operator is written as ci a = f% a & , where fi a 
is the spinon annihilation operator and e~ tei , the lower- 
ing operator of the 'angular momentum' Li which cor- 
responds to the charge quantum number. The enlarged 
Hilbert space is constrained by Li = J2 a f tafia — 1- The 
partition function is written as a path integral of e~ s ° 
over /, /*, 9 and h where 

S ° = J Q dT E {9t + iki ~ + 2U E( 9 ^ + ^) 



(tijfiafjat 



■c.c.) + ^(ji(l-x)-ihi)\.(2) 
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Here hi is the Lagrangian multiplier field imposing the 
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constraints Li — (J2 a flafia — 1) = 0, f-i is the chemical 
potential and x, the hole doping density from half filling. 
In the above action every term is quadratic in / and 9 ex- 
cept the hopping term. In Ref.[(| the hopping term was 
decomposed as a/3 «< a>P + (3<a> — < a >< (3 > 
with a = J2a fiafjv an d P = e l( - e '~ 6j \ Here we de- 
vise a more systematic way of decomposing the hopping 
term, so that the saddle point coincides with the previ- 
ous mean-field theory. The merit of the present method 
is that it enables us to expand about the saddle point to 
obtain a low energy effective theory. To decompose the 
hopping term we use the identity 



e - v , = / dr)i j dr) ij drijidr) ji 
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= Ea/it/jv and /% = 

g.^.-vj-Ay)^ This is a bit different from the conven- 
tional Hubbard-Stratonovich decomposition of a com- 
plete square and we encounter two problems. First, we 
need to introduce two independent complex fields r\ij 
and r]ji to decompose the hopping term and its com- 
plex conjugate, so that for a given configuration of r)ij 
and rjji, the action is not real. Second, it appears that 
the same rjij field decouples the fermion and the rotor 
term and it is not clear how we can recover the Florens- 
Georges's mean-field theory where < a >^< (3 >. Both 
difficulties are resolved as follows. We change the vari- 
ables of integration by rjij = \xtj \e Wil+llyai i +aij " > and 

Vji — \Xij ;\e~ Wl3+l<,a ^~ atJ \ Here \xij\ parameterizes the 
geometric mean of \r\ij | and \r/ji\, and my, the relative 
amplitudes, at and 8y are respectively the symmetric 
and antisymmetric combinations of the phases of rjij and 
rjji. It is emphasized that to !: j and a t are necessary in 
the integration because rjij and rjji are independent com- 
plex variables (not complex conjugate to each other). We 
follow Florens and Georges to replace e ldi by the boson 
variable Xi with a constraint |JQ| = 1 which is imposed 
by a Lagrangian multiplier A^. The partition function 
now takes the form 

Z = J D\xij\Dw ij Da ij Da±DfDf*DXDX*DhD\ 
[n <ij> , T (4| X y| 3 )]e- s , (4) 

where 

S= Jt 2 V \Uj\\x l3 \ 2 cosh2wij 
J ° 1 <id> 

+ E -x) -ihi +y,/* a (d T +ihi - /u)/ i(7 
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J2i( id r + hi)X*} [(-id T + hi)Xi] + I J2 W 
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3 — Wij+iatj +iAij j£* j£ 



Note that a.y and hi correspond to the gauge field associ- 
ated with the U(l) gauge transformation : — * e %Vi ji a , 
Xi — > e^X;, Oij — > a y - + (i/3, - </Jj) and /i, — > hi - d T ifi. 
Here we choose to couple the external EM field (Aij) to 
the boson. Alternatively can be coupled to spinon 
by a shift of the gauge field a,j — > a„ — Ay . 

The action © remains real for the fluctuations of \xij\ 
and a^. It becomes complex for the fluctuations of Wy, 
at, Xi and /ij. For the latter set of variables we lift the 
contour of integration into complex plane via analytic 
continuation to find a saddle point in the imaginary axis. 
Allowing for fluctuations along the deformed contour we 
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write \Xij 



m\ 



$\Xij\ 

+ 



Sai 



Wi 



IWi 



Sat, A 



J V i a ij = l ®tj + <* a ij> ^» = + ^ aild = + Shi, 

where quantities with tildes are real and represent the 
saddle point value. The appearance of imaginary saddle 
point is familiar in the treatment of the constraint field in 
the usual slave-boson models. In particular, fluctuations 
in the Sat and dwij directions are stable and massive. 
Together with 5|xij'|) these massive modes are negligible 
for energy scale less than t or U. The low energy effective 
Lagrangian is reduced to fermions and X bosons coupled 
to the compact U(l) gauge fields and hi, together 
with a constraint field Xi . If we choose, we can integrate 
over the Xi field to restore the 8 field resulting in the 
effective Lagrangian 

L = E fja(9r - ia T j - hi - fj,)f j<7 

3-P 

±- E(^ - < + Ik) 2 - E \Ux*e ia ?f; +a j ja 



2U 



,-i(e j+a -e:j-af-Af) 
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where a refers to the direction of neighboring sites con- 
nected by the hoppings. a" and aj are respectively the 
spatial and temporal gauge fields coming from 5a,j+ ai j 

and Shj in JSJl. The saddle point {x{j, xfj, hi, Xi) is 
determined from the extremum condition of the free en- 
ergy and is identical to the mean-field theory of Flo- 
rens and Georges 6], where x{j = \xij\e ai i~ lly ' w, - 3 ~ a *^ and 

x ^ = \xi \e~ a ^ +l( - Wl3+aij) . It is interesting to note that 
the auxiliary fields associated with the amplitude (lUy) 
and phase (at) fluctuations in the original integration 
@ switched their roles to make the free energy real at 
the saddle point and resolve both difficulties mentioned 
earlier. Eq. © is the first main result of this paper. 
, Here we remark on the connection between the 
slave-rotor theoryQ and the conventional slave-boson 
theory 0. In the slave-rotor theory the boson is rel- 
ativistic while the conventional slave-boson fields are 
non-relativistic. The relativistic X boson contains 
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particle and anti-particle. In the real-time canonical 
quantization the boson operator is written X(x,t) ~ 

/ % (b k e- l< - Ekt - kx 1 +d\e^ Ekt - kx ^, where the b k (d k ) 

is the (anti) particle annihilation operator and E], is the 
boson energy dispersion. b k (dk) carries positive (nega- 
tive) charges for both external (Aij) and internal (a^) 
gauge fields. Thus we identify b\ as creating holon and 
d\ as doublon. In the insulating phase the bosons are 
gapped, and holon and doublon are bound. In the metal- 
lic phase they becomes gapless with unbound holon and 
doublon. 

We apply the above formalism to find mean-field states 
in the triangular and honeycomb lattices at zero tempera- 
ture. We focus on the half filling where h — 0. For the tri- 
angular lattice we consider the flux phase : x(j — X^ e j 

and xfj = X Xe ^ j i where the phases are arranged so that 
the spinon (boson) acquires a flux when it moves around 
the smallest triangular plaquette. Here we choose t > 0, 
which can be always made by a particle-hole transforma- 
tion at half-filling. We have checked that the staggered 
flux states, where the fluxes take on alternating signs, 
are alway unstable relative to the zero flux state in the 
parameter range studied. For the Heisenberg model on 
the triangular lattice the fully gapped state with 7r-flux 
through the unit cell (henceforth, '7r-fluxphase') is known 
to be stable among paramagnetic states 8] . This is a state 
with flux 7r/2 through each triangle; it breaks time re- 
versal symmetry and is closely related to the chiral spin 
stated . It is of interest to see whether a gapless state 
without flux (henceforth, 'uniform phase') can be stabi- 
lized by charge fluctuations near the Mott transition. In 
our mean field theory the metallic phase is characterized 
by a nonvanishing amplitude of bose condensation which 
is given by Z — < X*(k — 0, r)X(k — 0, r) >, where 
X{k,r) is the boson field in the momentum space and 
N, the number of unit cells. For the isotropic triangular 
lattice a first order transition to the 7r-flux phase is found 
to occur at U c /t w 2.7 as is shown in Fig. 1. From finite 
size scaling we found that the uniform phase is stable 
over 2.6 < U/t < 2.73 as is shown in the inset of Fig. 
1. This insulating phase corresponds to the spin liquid 
state with a Fermi surface of neutral spinon. In order to 
see how the spin liquid phase is affected by anisotropy in 
the hopping integral we consider t ^ t with t , the hop- 
ping integral in one anisotropic direction. As \t — t\/t 
increases the region of the spin liquid phase shrinks and 
closes at \t — t\/t « 0.1 . In comparison with the numer- 
ical result |5j there are quantitative discrepancies. For 
example, Imada et al.[j| found U c /t w 5 for t = t and 
the spin liquid state is found to be gapless and therefore 
not the chiral state already for t jt = 0.5. The mean field 
phase diagram is not expected to be quantitative and it is 
apparent that the stability of the 7r-flux phase has been 
overestimated. Rather, the usefulness of our approach 
lies in obtaining a description of the insulating state and 
we now proceed to analyse the uniform phase. 



In the insulating phase the rotor has a gap and can 
be integrated out. This generates the Maxwellian kinetic 
term, C = -^(d x a) 2 with the gauge coupling, g 2 ~ ^- 

with A, the charge gap and n ~ tx? , the bare phase stiff- 
ness of boson [ljj . The low energy theory is reduced to the 
compact U(l) gauge theory coupled with spinons with 
Fermi surface. The compact U(l) pure gauge theory is al- 
ways confining but whether deconfinement is possible 
in the presence of matter field is an open question. Her- 
but et al. has argued that the theory is always confining 
in the presence of Fermi surface |l2^ or nodal fermions |l3j. 
Their conclusions actually apply to an approximate ef- 
fective action for the gauge field obtained by integrating 
out the fermions. Integrating out the fermions is sus- 
pect in the presence of gapless fermions. Indeed recently 
Hermele et al. |l ll proved that if the spin index is gener- 
alized to N flavors, the problem of 2N 2-component Dirac 
fermions coupled to complex U(l) gauge fields is decon- 
fined for sufficiently large N, thus providing a counter 
example to Ref. 13]. The spinon Fermi surface contains 
even more low energy excitations. In our view it is likely 
that a deconfinement state is possible and we will proceed 
with the assumption. 

There is extensive work on the coupling of the Fermi 
sea with non-compact gauge fields. A controlled per- 
turbation in l/N expansion is possiblefTEl Il6l Il7|. For 
small q and u, a quantum Boltzmann description is 
possible, albeit with singular Landau parameters which 
lead to anomalous power law behavior of many physi- 
cal quantities [Is|. For q = 2kp, there are singularities 
in the vertex function and it remains open whether the 
spin susceptibility is divergent or not|l9|. In contrast to 
the square lattice, there is no Fermi surface nesting for 
the half filled triangular lattice and antiferromagnetic or- 
der does not emerge naturally in our mean-field theory. 
Thus we believe that the spinon Fermi surface is stable 
and is the likely candidate for the spin liquid state found 
numerically |5j as well as experimentally [l|. Our conclu- 
sion is in agreement with that of Motrunich|3j . Of course, 
there remain the possibility of pairing instability due to 
residual interaction in any Fermi liquid, but that may set 
in at a very low energy scale. 

What are the physical implications of the spinon Fermi 
surface? Low lying gauge fluctuations produce a T 2 / 3 
contribution to the specific heatQ- The spinons carry 
spin and entropy. The spin conductivity (which is diffi- 
cult to measure) is proportional to the spinon transport 
time due to scattering by gauge fluctuations which goes 
as T~ 4//3 ^(|. On the other hand, the thermal conductiv- 
ity is proportional to the energy relaxation time which 
goes as T~ 2 / 3 because small q gauge fluctuations are ef- 
fective in relaxing the nonequilibrium energy distribu- 
tion. Using the quantum Boltzmann description^^] we 
predict that n/T ~ T~ 2 / 3 in the clean limit. The diver- 
gent behavior will be cut off by impurity scattering. Of 
course, being an insulator the charge conductivity is zero. 
This gross violation of the usual Wiedermann- Franz law 
together with unusual power laws are unique signatures 
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of this particular spin liquid state. In view of the un- 
usual behavior of K- (BEDT - TTF) 2 Cu 2 (CN) 3 it will 
clearly be of great interest to perform specific heat and 
thermal conductivity measurements in these materials. 

We have also applied the mean field theory to the hon- 
eycomb lattice. There is no obvious experimental candi- 
date but the interest in the honeycomb lattice is mainly 
theoretical. For nearest neighbor hopping, the fermion 
dispersion is characterized by 2 inequivalent Dirac cones 
at the Brillouin zone corner. On the insulating side of 
the Mott transition, the rotor is again gapped and the 
problem reduces to 2N Dirac cones coupled to a compact 
U(l) gauge field with N = 2. As mentioned earlier, there 
is hope that the theory is deconfined0|- Furthermore, 
the linear density of states greatly reduces the spin sus- 
ceptibility at a wavevector connecting the modes, and 
the AF ordered state does not emerge naturally in our 
mean field theory. It is amusing to observe that the en- 
ergy per bond of the naive AF state (— J/4) is identical 
to that of the dimer trial wavefunction. In this sense the 
honeycomb lattice is between the one dimensional chain 
and the square lattice. Our hope is that the AF state 
is not the ground state near the Mott transition, even 
though we can not address the competition directly in 
our mean field theory. Then the remaining competition 



is between the dimer state (which breaks rotational sym- 
metry) and the nodal spin liquid. Although not shown 
here we find for the nearest neighbor Hubbard model the 
nodal spin liquid state is stable over the narrow range 
1.68 < U/t < 1.74 near the Mott transition. Here we add 
further hopping terms t and t as shown in Fig. 2. It is 
straightforward to show that t does not affect the phase 
diagram or the energy as long as it is small compared to 
t. For negative t /t, we found that the spin liquid state 
is stable over a wider range of U/t, e.g., 2.3 < U/t < 2.6 
for t /t — —0.4 as is shown in Fig. 2. This region is the 
spin liquid phase with nodal spinons. Unusual power law 
for the spin response function correspon ding to AF fluc- 
tuations has been proposed for the state [2(|. Unlike the 
triangular lattice, the Hubbard model on the honeycomb 
lattice with t and t can be studied by quantum Monte 
Carlo without the Fermi sign problem so that a definite 
answer should be attainable|2ll Y22\ . 
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FIG. 1: The amplitudes of hopping order parameters 
of spinon (x^) and boson {x X ), the bose condensation 
amplitude(Z) and the uniform flux ($ in unit of n) per 
unit cell in the isotropic triangular lattice at half filling and 
T = OK. 50 x 50 lattice is used for the figures in the main 
pannel. The inset shows the size dependence of Z near the 
Mott transition. The arrow indicates the direction of increas- 
ing size from N = 20 x 20 to 80 x 80. The thick line denotes 
Z in the thermodyamics limit obtained from the finite size 
scaling. 
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FIG. 2: The amplitudes of the nearest neighbor hopping or- 
der parameters of spinon and boson^^), and the bose 
condensation amplitude(i?) in the honeycomb lattice with 
t'/t = and t" /t = -0.4 at half-filling and T = OK. The 
subscripts 1 and 2 refers to two inequivalent bonds in the 
dimerized phase. 



